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Abstract 

In this article we prove convergence of adaptive finite element methods for second order elliptic 
eigenvalue problems. We consider Lagrange finite elements of any degree and prove convergence 
for simple as well as multiple eigenvalues under a minimal refinement of marked elements, for all 
reasonable marking strategies, and starting from any initial triangulation. 

1 Introduction and main result 

In many practical applications it is of interest to find or approximate the eigenvalues and eigenfunctions of 
elliptic problems. Finite element approximations for these problems have been widely used and analyzed 
under a general framework. Optimal a priori error estimates for the eigenvalues and eigenfunctions have 
been obtained (see [B091, B089, RT83, SF73] and the references therein). 

Adaptive finite clement methods arc an effective tool for making an efficient use of the computational 
resources; for certain problems, it is even indispensable to their numerical resolvability. A quite popular, 
natural adaptive version of classical finite clement methods consists of the loop 

Solve — * Estimate — > Mark — > Refine, 

that is: solve for the finite element solution on the current grid, compute the a posteriori error estimator, 
mark with its help elements to be subdivided, and refine the current grid into a new, finer one. The 
ultimate goal of adaptive methods is to equidistribute the error and the computational effort obtaining 
a sequence of meshes with optimal complexity. Historically, the first step to prove optimality has always 
been to understand convergence of adaptive methods. A general result of convergence for linear problems 
has been obtained in [MSV07], where very general conditions on the linear problems and the adaptive 
methods that guarantee convergence are stated. Optimality for adaptive methods using Dorfler's [D6r96] 
marking strategy has been proved in [CKNS07, Stc07] for linear problems. 

The goal of this article is to analyze the convergence of adaptive finite element methods for the 
eigenvalue problem consisting in finding A £ R, and u ^ such that 

-V • (AVu) = XBu in 0, u = on 90, 

under general assumptions on A, B and f2 that we state precisely in Section 2.1. 

As we mentioned before, adaptive methods are based on a posteriori error estimators, that are 
computable quantities depending on the discrete solution and data, and indicate a distribution of the 
error. A posteriori error estimators for eigenvalue problems have been constructed by using different 
approaches in [Ver96, Ver94, DPR03, LarOO], they have been developed for A = I and B = 1, but 
the same proofs can be carried over to the general case considered here; see [GG07] and Section 2.3. 
An important aspect to be mentioned here is that the upper bound holds for sufficiently fine meshes. 
However, our proof will not rely on this bound, allowing us to prove convergence from any initial mesh. 
The first result (and only up to now) about convergence of adaptive finite elements for eigenvalue problems 
has been presented in [GG07]. 

The following is the main result of this article. 
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Theorem 1.1 (Main Result). Let Xk and u% be the discrete eigenvalues and eigenf unctions obtained with 
the adaptive algorithm stated in Section 2.4 below. Then there exists an eigenvalue X of the continuous 
problem such that 

lim Afe = A and lim dist H i m)(uk, M(X)) = 0, 

k— >oo k^oo 

where M(A) denotes the set of all eigenfunctions of the continuous problem corresponding to the eigen- 
value A. 

Remark 1. Before proceeding with the details of the statement and the proof of this result, we make 
some remarks: 

• An important difference with previous works is that we do not require the initial mesh % to be 
fine enough. Any initial mesh that captures the discontinuities of A will guarantee convergence. 

• The result holds for any of the popular marking strategies, not only Dorfler's [D6r96]. The only 
assumption is that non-marked elements have error estimators smaller than marked ones, see 
condition (2.17) in Section 2.4 below. 

• The marking is done according to the residual type a posteriori error estimators presented in 
Section 2.3. Even though there are some oscillation terms in the lower bound, we do not require 
any marking due to these terms. We only need to mark according to the error estimators, which 
is what is usually done in practice. 

• The result holds with a minimal refinement of marked elements, one bisection suffices. We do not 
require the enforcement of the so-called interior node property. 

• We are assuming that each of the discrete eigenvalues Xk is the j-th eigenvalue of the corresponding 
discrete problem. The result, as stated above, only guarantees that A^ converges to one eigenvalue 
A of the continuous problem. We can be sure that we approximate the j-th eigenvalue of the 
continuous problem under any of the following assumptions: 

— No eigenfunction is equal to a polynomial of degree < t on an open region of f2, where i 
denotes the polynomial degree of the finite clement functions being used. This is a Non- 
Degeneracy Assumption, and it holds for a large class of problems; see Assumption 1 and 
following discussion. 

— The meshsize of the initial triangulation is small enough. This assumption goes against the 
spirit of adaptivity and a posteriori analysis, since we cannot quantify what small enough 
means. But we state it for completeness, because in some (nonlinear) problems there may be 
no way to overcome this. 

• The proof follows similar ideas to those of [MSV07], with some modifications due to the different 
nature of the problem. It consists in proving the following steps: 

— The full sequence of discrete eigenvalues converges to a number Aoo and a subsequence of the 
discrete eigenfunctions converges to some function Uoo- 

— The global a posteriori error estimator converges to zero (for the subsequence). 

— The pair (Aoo, Moo) is an eigenpair of the continuous problem. Due to a lack of a sharp upper 
bound (it only holds for sufficiently fine meshes) it is necessary to introduce a new argument to 
prove this (see Theorem 5.1). This new argument is perhaps the main difference with respect 
to [MSV07], and we believe that the idea can be useful for many other nonlinear problems. 

The rest of the article is organized as follows. In Section 2 we state precisely the problem that we 
study, describe the approximants, mention some already known results about a priori and a posteriori 
estimation, and state the adaptive loop. In Section 3 we prove that the sequence {(Xk, Uk)}keN °f 
solutions to the discrete problems contains a subsequence that converges to a limiting pair (Aoo, ^oo)- In 
Section 4 we prove that the global a posteriori error estimator tends to zero; which is instrumental to 
conclude in Section 5 that (Aoo,Uoo) is an eigenpair of the continuous problem. Finally, in Section 6 we 
state and prove the main result and discuss its implications. 
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2 Problem statement and adaptive algorithm 



This section is subdivided in four parts. In Section 2.1 we state precisely the continuous problem that 
we study and mention some of its properties. In Section 2.2 wc state the discrete problems that we 
consider as approximants to the continuous one, mention some of its properties and state the a priori 
error estimates. In Section 2.3 we define the a posteriori error estimators that we use, state the upper 
bound and prove the discrete local lower bound that we will use in our convergence proof. Finally, in 
Section 2.4 we state the adaptive algorithm together with the assumptions on each of its blocks. 

2.1 Setting 

Let f2 C M. d be a bounded open set with a Lipschitz boundary. In particular, we suppose that f2 is a 
polygonal domain if d = 2 and a polyhedral domain if d = 3. Let a,b : Hq(CI) x Hq(SY) — > R be the 
bilinear forms defined by 



a(u,v) := I AVu-Vv, 

and 



b(u, v) := / Buv 



Q 

where A is a piecewise W 1,oa (£l) symmetric-matrix- valued function which is uniformly positive definite, 
i.e., there exist constants a\, a 2 > such that 

ailCI 2 < A(x)£ ■ £ < a 2 \£,\ 2 , V£eR d , V x E O, 

and B is a scalar function such that 

h < B(x) <b 2 , VxeSJ, 

for some constants b%, 62 > 0. 

We also define the induced norms by these bilinear forms as 

\\v\\ a :=a{v,vfl 2 , vEH^n), and \\v\\ b := &(^) 1/2 , v E L 2 (n). 

By the assumptions on A and £>, ||-|| a ~ || • ||ffi(o) and ||-|| fc ~ || • i.e., there exist positive constants 
Ci, C2, C3, C4 such that 

Ci\\v\\ Him < \\v\\ a < c 2 \\v\\ H i {n} , VvEH^n), 

and 

c 3 \\v\\a< IHI 6 <C4|H|n, V«eL 2 (Sl). 

Where, hereafter, if A C ft, \\ ■ \\a will denote the I/ 2 (A)-norm. 
We consider the following 

Continuous eigenvalue problem. Find Ael and u E Hq(V£) satisfying 

a(u,v) = Xb(u,v), V v E H^(n), , 
\\u\\ b = l- (2A) 

It is well known [B091] that under our assumptions on A and B problem (2.1) has a countable 
sequence of eigenvalues 

< Ai < A 2 < A 3 < . . . / 00 

and corresponding eigenfunctions 

U\,U2, u 3 , . . . 

which can be assumed to satisfy 

1 i=j, 



b(ui,Uj) = 8 l3 



i+h 
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where in the sequence {Aj}j £ n, the Xj are repeated according to geometric multiplicity. 

a{u, tx) 

Also, the eigenvalues can be characterized as extrema of the Rayleigh quotient lZ(u) = — — ■ — ^, by 



the following relationships. 
• Minimum principle: 



b(u, u) 



Ai = min 7Z(u) = TZtui), 

A, = min lZ(u) = 1Z(uj), j = 2, 3, 



• Minimum-maximum principle: 



Xj = min maxTZ(u) = max 

VjdHoityuZVj u6span{«!,. 
dim Vj —j 



K{u), 3 = 1,2,. 



For each fixed eigenvalue A of (2.1) we define 

M(A) := {u G H^(Q) : u satisfies (2.1)}, 

and notice that if A is simple, then M(A) contains two functions, whereas if A is not simple, it consists 
of a sphere in the subspace generated by the eigenfunctions. 



2.2 Discrete problem 

In order to define the discrete approximations we will consider conforming triangulations T of the domain 
n, that is, partitions of O into d-simplices such that if two elements intersect, they do so at a full 
vertex/edge/face of both elements. For any triangulation T, S will denote the set of interior sides, where 
by side we mean an edge if d = 2 and a face if d = 3. And Kr will denote the regularity of T, defined as 

diam(T) 
Kt '■= max , 

where diam(T) is the length of the longest edge of T, and pt is the radius of the largest ball contained 

in it. It is also useful to define the meshsizc hj- '■= max /it, where hx '■= |T| 1//d . 

Ter 

Let I £ N be fixed, and let Vt be the finite clement space consisting of continuous functions vanishing 
on dfl which are polynomials of degree < £ in each element of T, i.c, 

V r := {v £ H^(Q) : v\ T e V e (T), VTe T}. 

Obviously, Vr C Hq(&1) and if % is a refinement of T, then Vr C Vr»- 

We consider the approximation of the continuous eigenvalue problem (2.1) with the following 



Discrete eigenvalue problem. Find At £ K and i*t G Vt such that 

( a(u T ,v) = X r b{u T ,v), VugVt, ^ 2 x 

I ||«r|| 6 = L 

For this discrete problem, similar results to those of the continuous problem hold [B091]. More 
precisely, problem (2.2) has a finite sequence of eigenvalues 

< Ai,t < A 2: t < A 3; t < ■ ■ • < Xn t ,t, 
where Nq- := dimVT, and corresponding eigenfunctions 

Ul.T, U2,T, «3,T, • ■ ■ , UN T .T, 

which can be assumed to satisfy 

6(u^t,Wj,t) = hj- 
Moreover, the following extremal characterizations also hold: 
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• Minimum principle: 

Ai.r = min lZ(u) = 7Z(ui.r), 

XjT = niin 1Z(u) = lZ(uj,r), j = 2, 3, . . . , Nq-. 

a(u,Ui t)—0 

i=i,.:.j~i 

• Minimum-maximum principle: 

XjT = min max lZ(u) = max 7t(u), j = 1, 2, . . . , JV7-. 

dim Vj,r=J 

It follows from the minimum-maximum principles that 

A, • A, 7. j = l,2,...,N T . 
and it also follows that if 71 is any refinement of T then 

A,./. • A,. 7 . ./ 1.2 V/. 

For a given eigenvalue A we define a notion of minimal error of approximation of its eigenfunctions 



by 

For j = 1,2,..., Nt, there holds that 



er(A) := sup inf ||u - x\\ a ■ 

uGM(A) XGV r 



Kr-^ <4(A,), 



where, from now on, whenever we write A < £> we mean that A < Ci? with a constant C that may 
depend on A, 23, the domain f2 and the regularity nq- of T, but not on other properties of T such as 
element size or uniformity. 

If {7fe}fc g N is any sequence of triangulations such that sup Kq- k < 00, and hq- k — > as k — * 00, then 

fc£N 

er fe (Aj) — > 0, as fc — > exo, 

and therefore, 

Aj,r fc — *■ Aj, as fc — > 00. (2-3) 

This holds for any j 6 N and it is a consequence of standard interpolation estimates and the fact that 
M (Aj) is bounded and contained in a finite dimensional subspace of H^fl). 



2.3 A posteriori error estimators 

A posteriori estimates for eigenvalue problems have been studied by Larson [LarOO], Duran, Padra and 
Rodriguez [DPR03] , Giani and Graham [GG07] . In this section we present the residual type a posteriori 
estimates for eigenvalue problems, state without proof some already known properties and prove the 
discrete local lower bound that will be useful for our convergence proof. 

In order to define the estimators we assume that the triangulation T matches the discontinuities of 
A. More precisely, we assume that the discontinuities of A are aligned with the sides of T. Observe that 
in particular, A\t is Lipschitz continuous for all TeT. 

Definition 1 (Element residual and jump residual). For /.t E R and v G Vr we define the element 
residual R(n, v) by 

R(ji,v)\ T ■= -V- (AVv) -iiBv, (2.4) 
for all T £ T, and the jump residual J(v) by 

J(v)\s := (AS7v)\ Tl ■ ni + (AVv)\t 2 ■ na, (2.5) 

for every interior side S £ S, where Ti and Ti are the elements in T which share S and nl is the outward 
normal unit vector of Ti on S, for i = 1,2. We define J(v)\gn := 0. 
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Definition 2 (Local and global error estimator). For /j 6 1 and v G Vr we define the local error 
estimator ijq-(fJ,,v;T) by 

ifr((J-,v;T) 2 := h\ \R{\i,v)\ T + h T ||J(i;)||| T , 
for all TsT, and the global error estimator t7t(a*; v ) i s given by 

Vt(^,v) 2 := ^2 Vt{^,v;T) 2 . 
TeT 

Even though we will not need it for the convergence proof, we include the statement of the upper 
bound of the error in terms of the a posteriori error estimation, for the sake of completeness. 

Theorem 2.1 (Upper bound). Let j G N, and let uq- be an eigen] 'unction corresponding to the j-th 
eigenvalue At of the discrete problem (2.2), then, if Jit is small enough, there exists an eigenfunction u 
corresponding to the j-th eigenvalue A of the continuous problem (2.1) such that 

\\u - UT\\ a < i]t{^t,u t ). 

The proof of this theorem can be obtained following the steps given in [DPR03] , by extending Lemmas 
3.1 and 3.2 presented there for the model problem with A = I, and B = 1, to the general case that we 
consider here, using the following regularity result, and the a priori bound stated in Theorem 2.2 below. 

Lemma 1 (Regularity of the cigenfunctions). There exists r G (0,1] depending only on D, and A such 
that 

u G H 1+r (Q), 

for any eigenfunction u of the problem (2.1). 

Proof. This can be proved by observing that if u is an eigenfunction, then it is also a solution to a linear 
elliptic equation of second order with right-hand side in L 2 (fl). We know that r — 1 when A is constant 
or smooth and O is convex. The case in which f2 is non-convex has been studied in [BDLN92] and the 
case of A having a discontinuity across an interior interface in [Bab70]. For the general case, which we 
are considering here, see [Joc99, Theorem 3]. □ 

The following result is an a priori estimate relating the errors in the strong and weak norms associated 
to the problem, and it is the last slab in the chain necessary to prove Theorem 2.1. The case A = I and 
6 = 1 can be easily obtained from the results in [SF73], and in [RT83]. The general case was presented 
in [GG07, Gia07]. 

Theorem 2.2. Let the same assumptions of Theorem 2.1 hold. Then, ifhq- is small enough, there exists 
an eigenfunction u corresponding to the j-th eigenvalue A of the continuous problem (2.1) such that 

\\ u - u r\\ b < h r T ||u-ur|| . 

The next result, which we will need for our proof of convergence is the discrete local lower bound, 
whose proof follows that of the continuous lower bound in [DPR03], but in order to make this article 
more self-contained we will include it here. 

For S G S we define ujt(S) as the union of the two elements in T sharing S. For TeT, Mt(T) := 
{T 1 G T : V n T ^ 0} denotes the set of neighbors of T in T, and w T (T) := Ur'GJVV(r) T '■ Wc also 
define := 3 if d = 2 and := 6 if d = 3. This guarantees that after na bisections to an element, 
new nodes appear on each side and in the interior. Here we consider the newest-vertex bisection in two 
dimensions and the procedure of Kossaczky in three dimensions [SS05] . 

Theorem 2.3 (Discrete local lower bound). Let T G T and let T„ be the triangulation of fl which is 
obtained from T by bisecting n c i times each element of Mq-{T). Let At and uq- be a solution to the 
discrete problem (2.2). Let W be a subspace of iJg(fi) such that Vt, C W. If fx 6l and w G W satisfy 

J a(w,v) = fib(w,v), V»eW, 

I IH 6 = i, 



G 



then 

t)t(\t,ut;T) < \\V(w -u r )IL r(T) +h T \\/iw - Xtut\\ Ut{t) + h T \\R - R\\ ut[t) + h^ 2 \\j - j\\ dT , 

where, for every T' £ J\fq-(T), R\t> is the L 2 (T') -projection of R := R(\t,ut) onto Vg-i, and for every 
side S C dT, J\s is the L 2 (S)-projection of J := J(ut) onto Vt—i- 

Proof. Q~] We first analyze the element residual. Wc obviously have 

\\R\\ T <\\R\\ T + \\R-R\\ T . (2.6) 

Let XJ? denote the vertex of T„ which is interior to T. Let ipx be the continuous piecewise linear function 
over % such that ifT = 1 and ipr vanishes over all the others vertices of %. Then 

Pllr~ / R 2( Pt= f R(R<p T )= I R(Rtp T )+ [ (R-RjRyT- (2.7) 
Jt Jt Jt Jt 

If we define v := Ripr G Vr, C W, taking into account that v vanishes over dT, for the first integral 
in (2.7) we have 

I Rv= I (-V • {AVu T ) - X T Bu T )v 
Jt Jt 

AVut ■ Vv — / XtButv 
Jt 

= I AVu r • Vu - / X T Bu r v - / AVw • V« + / \x&wv 
Jt Jt Jt Jt 

= / AV(ut — w) ■ Vu + / B(fj.w — Xtut)v 
Jt Jt 

< ||V(ur — w)\\t \\^ v \\t + \\t lw ~ XtutWt \\ v \\t ■ 
For the second integral in (2.7) we have 

/ (r — R)R(fT < p^T|Lp-i?|L < p|Lp--R|L. 

Jt 

Therefore, taking into account that ||Vw|| T < \\v\\ T and ||k|[ t < P|| T wc can write 
PIIt ~ II V(ut - w)\\ T — \\R\\ T + \\/J,w - Xtut\\ t \\R\\ T + \\R\\ T \\R — R\\ T , 

and then 

hr p|| T < ||V(u T - w)\\ T + h T \\fxw - Xtut\\ t + h T \\~R~ - R\\ T . (2.8) 
Now, from (2.6) and (2.8) it follows that 

h T \\R\\ T £ ||V(tir - u>)|| T + h T \\nw - Xtut\\ t + h T \\R - ~R~\\ T . (2.9) 
The same bound holds replacing T by T", for all T 6 JV r (T). 

Secondly, we estimate the jump residual. Let S be a side of T and let T\ and T 2 denote the elements 
sharing S. Obviously, one of them is T itself. As before we proceed by bounding first the projection J 
of J, since 

\\J\\s <\\J\\ S + \\J -J\\ s - (2-10) 

Let Xg* denote the vertex of %, which is interior to 5". Let ips be the continuous piecewise linear function 
over % such that ipsfa™*) = 1 and ips vanishes over all the others vertices of %. Then 

s 



JL< / (J)Vs= / J(Jvs)= / J(J<Ps)+ / (J-J)J^s- (2.11) 
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Now, we extend J to u>t(S) as constant along the direction of one side of each T;, for i = 1,2, and 
still call this extention J. Observe that J is continuous on u>t(S) and J|t; G / P^_i(T,), for i = 1,2. 

Since i; := J</?s G Vr» C W and taking into account that u = on d(ojr(S)), for the first integral 
in (2.11) we have 

/ J v = V / MVur • nt = Y] / V ' («-4Vur) 
= ^2 (^J AVu T ■ Vv + J uV • (-4Vit r )^ 

= V ( / -4Vu r -Vv+ uV • (iVar) ) + / - / • Vv 

i=l,2 ' JT1VJT2 JT 1 UT 2 

= (yj AV{u T -w)-Vv + J (V • (.AVut) + 



ML 



= / AV(u T - w) -Vv + V ( / -Rv+ / B(fj,w - \ r u r )v ) 

< ||V(ut - ™)IL r( s) ll Vu L r (s) + II^Lx(s) IMl-r(S) + H/^ - A r^rlL r(s) imi wr(s) 
For the second integral in (2.11) we have 

(J-J)J^ S < P<Ps\\ s P-J\\ s < IPIIsp- J\\ s - 



Hence, taking into account that ||Vu||^ (s) < |ML r(s) , 



h^ 2 \\j\\ s we can write 



Ml,,^^ < \\J\\ rcn an d K /m 



< 



\J\\ 2 s ^\M^t-w)\I hs) h T 1/2 \\j\\ s ^ 

+ \\flW - \TUT\\ UT (S) h T 2 \\ J \\s 



\R\ 



wt(S) ,l T 



h 1 ' 2 Oil 



Kllslr J \\s 



and then 



h-r 2 \\J\\ S < ||V(ur - w)IL T (5) + ^ II^ILrfS) + ^ - A ritrL r(s) + /ly 2 II J - J |ls 



(2.12) 



Now, from (2.10) and (2.12) it follows 



h T 2 \\ J \\s ~ ll V ( u T - W)\\ UT ( S ) + h T \\R\L T (S) + h T \\V W - X TU T \L t{S ) + 



1/2 



Adding the last equation over all 5 C dT, we obtain 

h r 2 \\J\\ar ~ ll v ("r - ™)IL r(T ) + II-RIL^t) + h T \\v w ~ x tu t 
The claim of this theorem follows by adding this last inequality and (2.9) 



+ h 1/2 



J -J 



J — J 



ST 



□ 



The next result is some kind of stability bound for the oscillation terms, which will be useful to obtain 
the bound of Corollary 1 which is what will be effectively used in our convergence proof. 

Lemma 2. Under the assumptions of Theorem 2.3 there holds 

h T \\R - R\\^ {T) + hlj, 2 J|| ST < h T {2 + \ r )\\u r \\ H i {ulr(T)) . 

Proof. Q~] We first consider the term corresponding to the element residual. 



\R — R\\ T = 



< 



-V • (AVu T ) - XtBut + V • {AVu T ) + XtBut 



-V • (AVu r ) + V • (AVu r ) + \\Xt{Bu t - Bu T )\ 



(2.13) 
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where, as before, the bar denotes the -L 2 (T)-projection onto Tt-\(T). 

Let A T = (Afj) denote the mean value of A = (Aij) over the element T, and note that V- (A t Vut) = 

V • (A t \7ut). Thus, for the first term in the right hand side of (2.13) we have 



-V • (AVu T ) + V • (A\7u T ) = V • ({A T - A)Vu r ) - V • {{A T - A)Vu T ) 



< 



V • (A T — A)) • Vu r - (V • (A T - A)) ■ Vu T 
(A T - A) : D 2 u T - (A T - A) : D 2 u T 

< \\(\7 ■ {A T ~ A)) ■ Vu T \\ T + \\{A T ~ A) :D 2 u T \ 

< UWwut) ||Vut|| t + U T - A\\ L ~ {T ) II^VL 



Since \\A T — A\\l^{t) S ^tWAWw 1 (t); an inverse inequality leads to 



-V • (AVu T ) + V • (AVu T ) < MHiv^(t) ||Vur|| T + Ml-^llwn (t) ||L> 2 mt|| t < ||Vur 
For the second term in the right hand side of (2.13) we have 

\\\t(But - But) || t < \\^tBut\\ t < At ||ur|| T , 

and therefore, 

\\R-R\\ T Z(l + *T)\\u T \\m(T)- 
The same estimation holds for all elements in Afq-(T), and consequently, 



h T \\R-R\ 



wt(T) 



< h T (l + \ T )\\uT\\m(u> T {T))- 



(2.14) 



HH Next, we analyze the jump residual. Let S be a side of T and let T\ and T2 denote the elements 
sharing S. Again, if the bar denotes the L 2 (S')-projection onto Ve-i(S), it follows that 



\J-J\\s = 



(AVu T )\T t ■ nj - ^ (.AVur)!^ 



Using that (^4 Ti Vut)|t; • "1 = (^4 Ti Vur)!^ ■ Hi we have 



J — J 



£ ((^l-A T 0Vu T )| Ti . -nt- 2 ((^-Ari)V«r)k 



i=l,2 



i=l,2 



< 



Therefore, 



J] ((^-A T 0Vu r )| Ti -n? 
=1,2 

^ Z ||((«^-^ Ti )V«r)|T 4 -m|| s 

i=l,2 

< E -^ T4 lli-(S)HV«rh|| s 

i=l,2 

< E Ml^llv^(T 4 )^ 1/2 ||ur||j?iCTi) 

i=l,2 

< h^ 2 \\u T \\ H i (uT(s)) . 

2 y 2 11 1 1 

h T \\J ~ J\\ dT <h T \\uT-\\ H i^ r ^ T) y 



Adding (2.14) and (2.15) we obtain the claim of this lemma. 

As an immediate consequence of Theorem 2.3 and Lemma 2 the following result holds. 
Corollary 1 (Lower bound). Under the assumptions of Theorem 2.3 there holds 

Vt(Xt,U T ;T) < \\V(W-U T )\\ UT ( T) + /l T ||/HL-r(T) + h T{l + \T)\\u r \\ H i{uT(T)y 



(2.15) 
□ 
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2.4 Adaptive loop 

Our goal is to use an adaptive method to approximate the j-th eigenvalue and one of its eigcnfunctions, 
for some fixed j G N. From now on, we thus keep j G N fixed, and let A denote the j-th eigenvalue of 
(2.1) and u an cigenfunction in M(A). 

The algorithm for approximating A and M(X) is an iteration of the following main steps: 

(1) (Afc,tt fc ) :=SOLVE(V fc ). 

(2) { Vk (T)} TeTk := ESTIMATE^, Ufc ,T fe ). 

(3) M k := MARK({ % (r)} TeTfc ,r fc ). 

(4) T fe+ i := REFINE(T fe ,X fe ), increment k. 

This is the same loop considered in [MSV07] , the difference lies in the building blocks which we now 
describe in detail. 

If T k is a conforming triangulation of il, the module SOLVE takes the space Y k := Vr* as input 
argument and outputs the j-th eigenvalue of the discrete problem (2.2) with T = T k , i.e., X k := Xj,T k i 
and a corresponding cigenfunction u k G Y k - Therefore, X k and u k satisfy 

( a(uk,Vk) = X k b(v,k,Vk), VufeGVfc, (2 16) 

I IKII 6 = 1- 

Given T k and the corresponding outputs X k and u k of SOLVE, the module ESTIMATE computes and 
outputs the a posteriori error estimators {rj k (T)}T£T k , where 

Vk(T) := r)T k (X k ,u k ;T). 

Based upon the a posteriori error indicators {r] k (T)}T^T k , the module MARK collects elements of T k 
in M. k . In order to simplify the presentation, the only requirement that we make on the module MARK 
is that the set of marked elements Ai k contains at least one element of T k holding the largest value of 
estimator. That is, there exists one element T" lax G M. k such that 

%(^ nax ) = max % (T). 

Whenever a marking strategy satisfies this assumption, we call it reasonable, since this is what practi- 
tioners do in order to maximize the error reduction with a minimum effort. The most commonly used 
marking strategies, e.g., Maximum strategy and Equidistribution strategy, fulfill this condition, which is 
sufficient to guarantee that 

/ \ 1/2 

TeT k \M k Vk(T) < Vk(Mk) ■= f J2 ^{T) 2 ] . (2.17) 

This is slightly weaker, and is what we will use in our proof. The original Dorfler's strategy also 
guarantees (2.17). 

The refinement procedure REFINE takes the triangulation T k and the subset M k C T k as input 
arguments. We require that all elements of A4 k arc refined (at least once), and that a new conforming 
triangulation T k+ \ of Q, which is a refinement of 7~ k , is returned as output. 

In this way, starting with an initial conforming triangulation % of f2 and iterating the steps (1),(2),(3) 
and (4) of this algorithm, we obtain a sequence of successive conforming refinements of Tq called T\ , T%, . . . 
and the corresponding outputs (X k , u k ), {i] k (T)}T£T k , -Mk of the modules SOLVE, ESTIMATE and MARK, 
respectively. 

For simplicity, we consider for the module REFINE, the concrete choice of the newest vertex bisection 
procedure in two dimensions and the bisection procedure of Kossaczky in three dimensions [SS05] . Both 
these procedures refine the marked elements and some additional ones in order to keep conformity, and 
they also guarantee that 

K := sup Kq- k < CO, 
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i.e., {Tk}keN is a sequence shape regular of triangulations of f2. It is worth mentioning that we do not as- 
sume REFINE to enforce the so-called interior node property, and convergence is guaranteed nevertheless, 
this is an important difference with respect to [GG07]. 

Regarding the module MARK, we stress that the marking is done only according to the error es- 
timators; no marking due to oscillation is necessary, this is another important difference with respect 
to [GG07], where the set of marked elements has to be enlarged so that Dorfler's criterion is satisfied not 
only by the data oscillation terms, but also by the oscillation of the current solution Uk- 

3 Convergence to a limiting pair 

In this section we will prove that the sequence of discrete eigenpairs {(Xk, Ufc)}fcgN obtained by SOLVE 
throughout the adaptive loop of Section 2.4 has the following property: Xk converges to some Aoo G M 
and there exists a subsequence {ufc m } me N of {iifc}fceN converging in i? 1 (fJ) to a function itoo- 

Let us define the limiting space as Voo := UVfc , and note that Voo is a closed subspacc of H^fl) 
and therefore, it is itself a Hilbcrt space with the inner product inherited from Hq(£1). 

Since Tk+\ is always a refinement of 7^, by the Minimum-Maximum principle {Afc}^^ is a decreasing 
sequence bounded below by A. Therefore, there exists Aoo > such that 

Afc \ Aqo. 

From (2.16) it follows that 

IWIo = o.{u k ,u k ) = X k b(u k ,u k ) = Xk \\uk\\l = X k -> Aoo, (3.1) 

and therefore, that {uk}keN is a bounded sequence in Voo. Then, there exists a subsequence {uk m } m eN 
weakly convergent in Voo to a function G Voo, so 

w fcm ^Uoo in H^ty. (3.2) 

Using Rellich's theorem we can extract a subsequence of the last one, which we still denote {wfe m }meN i 
such that 

Uk m — > "oo in L 2 (Q). (3.3) 

If ko G No and k m > ko, for all Vk G Vfc we have that a(uk m ,Vk ) = Xk m b(uk m ,Vk ),, and when m 
tends to infinity, we obtain that a(uoo,Wfe ) = Aoo6(«oo, v k a )- Since ko G No and Vk G Vfc are arbitrary 
we have that 

a(woo, v) = Aoo^Woo, v )> VwGVoo- (3.4) 
On the other hand, since that j|ufe m || fc = 1, considering (3.3) we conclude that ||uoo!| fc = 1- Now, taking 
into account (3.4) we have that 

|| u oo|| a = Aqo I W-oo 1 1 & = Aoo- 

From (3.1) it follows that ||iifc m ||^ = Xk m — > Aqo, and therefore, ||w/c m || a — > ||woo|L- This, together 
with (3.2) yields 

Uk m — ► Uoo in Hq(Q,). 
Summarizing, we have proved the following 
Theorem 3.1. There exist Aqo G K and u M G Voo such that 

( a(uoo, v) = Aoo biuoo, v), VuG Voo, 

I ll«oo|| t = l. 

Moreover, Aoo = hm A& and there exists a subsequence {uk m }meN °/{ M fe}feGN such that 

k — >oo 

Uk m — > "in Hq(Q). 

Remark 2. It is important to notice that from any subsequence {(Afe m , Uk m )}meN °f {(^fcj w fc)}fceN ) 
we can extract another subsequence {(Afe m ,«i, m )}neN ) such that Uk m converges in i? 1 (fi) to some 
function Uqo G Voo that satisfies 

f a(Uoo,v) — Aoo 6(iioo, w), VueVoc, 
1 ||«oo||i, = 1- 
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4 Convergence of estimators 

In this section wc will prove that the global a posteriori estimator defined in Section 2.3 tends to zero. We 
will follow the same steps as in [MSV07] providing the proofs of the results that are problem dependent. 
Those geometrical results that arc consequences of the fact that we are only refining will be stated 
without proof, but with a precise reference to the result from [MSV07] being used. 

In order not to clutter the notation, we will still call {itfc}fceN to the subsequence {uk m }meN , and 
{Tk}keN to the sequence {^fc m }meN - Also, we will replace the subscript 7^ by k (e.g. N k {T) := Nr k (T) 
and u>k{T) := u>r k (T)), and whenever S is a subset of r/k(E.) 2 will denote the sum X)tg~ Vk(T) 2 . 

The main result of this section is the following 

Theorem 4.1 (Estimator's convergence). If {Tk}keN n denote the triangulations corresponding to the 
convergent subsequence of discrete eigenpairs from Theorem 3.1, then 

lim Vk (T k ) =0. 

k — >oo 

In order to prove this theorem we consider the following decomposition of T~k, which was first estab- 
lished in [MSV07] . 

Definition 3. Given the sequence {Tk}ke^ °f triangulations, for each k G No wc define the following 
(disjoint) subsets of Tf.. 

• T k ° := {T G % : T is refined at least n d times, for all T G A4(T)}; 

• T+ := {T e T fe : T is never refined, for all V G M k (T)}; 
. T fe *:=T fc \(T fc °UT+). 

We also define the three (overlapping) regions in O: 
. fi° := U w fc (T); 

TGT L ° 

. 0+ := U 

r e r+ 

. Sl* k ■= U w fc (T). 

TGT fc * 

We will prove that T)k(T k °), rj k (T k *) and rj k (T k + ) tend to zero as fc tends to infinity in Theorems 4.2, 
4.3 and 4.4. Since rj k {%) 2 = r/ k (T k ) 2 + ?] k (T k + )' 2 + r] k (T k *) 2 , Theorem 4.1 will follow from these results. 

Definition 4 (Meshsize function). We define h k G L°°(f2) as the piecewise constant function 

h k \ T ■= \T\ 1/d , VTeTfc. 

For almost every x G fl there holds that h k (x) is monotonically decreasing and bounded from below 
by 0. Therefore, 

hoo{x) := lim h k (x) 

k — >oo 

is well-defined for almost every x G SI and defines a function in L°°(0). Moreover, the following result 
holds [MSV07, Lemma 4.3 and Corollary 4.1]. 

Lemma 3. The sequence {/ifc}fc<=N converges to hoo uniformly, i.e., 

lim \\h k - fooo||z,°°(n) = 0, 

k — >oo 

and if xn° denotes the characteristic function of Q k then 

lim \\h k Xn°\\L°°(n) = 0. 
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This lemma is a consequence of the fact that the sequence of triangulations is obtained by refinement 

only, and that every time an element T G T k is refined into "Tfe+i, hk+i(x) < i^) 1 ^ h^ix) for almost every 
x G T. But it is otherwise independent of the marking strategy. The next result is also independent 
of the marking strategy, it is just a consequence of the fact that u k — > u^, the lower bound and the 
convergence of \\hkXa° h IU°°(n) to zero. 

Theorem 4.2 (Estimator's convergence: First part). If {Tk}ke^ denote the triangulations correspond- 
ing to the convergent subsequence of discrete eigenpairs from Theorem 3.1, then the contribution of T k ° 
to the estimator vanishes in the limit, i.e., 

km r lk (T k °) = 0. 

ft— i-OO 

Proof. Using Corollary 1 with W = Yoo, w — ^oo ; and [i = Aqo we have that 

Vk(^) 2 = J2 ^ T ) 2 
TeT" 

^ £ W V ( U k - u °°)\\l k (T) + h T l|AooMoo||^ (T) + h 2 T {l + Afc) 2 |K||^i (wfc(T)) 

TeT° 

< ||V(lifc - "oo)i|j 2 + H^fcX^n ( IIAooMooH^ + (1 + Afc) 2 ||w fc ||^i (n) ) 

Since X k — > Aoo in R and — ► Uoo in £T (fi), Lemma 3 implies the claim. □ 

The following lemma was proved as the first step of the proof of [MSV07, Proposition 4.2], it is also a 
consequence of the fact that the sequence of triangulations is obtained by refinement, without coarsening, 
and it is independent of the specific problem being considered. 

Lemma 4. If Qj£ is as in Definition 3, then 

km |^|=0. 

K — >00 

From Corollary 1, Lemma 4 and tke fact that u k — ► en H 1 we obtain 

Theorem 4.3 (Estimator's convergence: Second part). If {T k } ke jq denote the triangulations corre- 
sponding to the convergent subsequence of discrete eigenpairs from Theorem 3.1, then the contribution of 
T k * to the estimator vanishes in the limit, i.e., 

km m {T k *) = 0. 

k — >oo 

Proof. Let (A, u) be any eigenpair of (2.1). Then Corollary 1 with W = Hq(CI), w = u, and fj, = A implies 
that 

miV) 2 = £ vk(Tf 

< £ l|v(ti fc - u)\\l k(T) + h% \\\u\\l k(T) + 4(i + ^) 2 \\u k \\ 2 m{uJk{T)) 

TeT k * 

< ||vk - u)\\ 2 ui + x 2 \\u\\ 2 a , + (i + A^HIIl^) 

< ||V(« fc - «oo)|| 2 2 + ||V(uoo - u)\\ 2 n + A 2 || u || 2 , 

+ (1 + \ k ) 2 \\u k - Uoollffifn) + (1 + Afc) 2 ||woo|||-i(n*)- 

Taking into account that X k — > A^ in R, u k — > t/qo in H l (Vl) and Lemma 4, the claim follows. □ 
In order to prove that the estimator contribution from T k + vanishes in the limit, we make the following 
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Definition 5. Let T + be the set of elements that are never refined, i.e., 

k>0 m>k 

and let the set Q + be defined as 

Q+:= |J T. 

TeT+ 

It is interesting to observe at this point that 
Lemma 5. The set f2 + is empty if and only if lim ||/ifc||.L°°(n) = 0. 

k— >oo 

Proof. If il + is empty, then fijj" and f2£ are empty for all k G No, and ||/ifc||i°°(f2) — ||^*]|z,«>(Q ) which 
tends to zero by Lemma 3. Conversely, if liirifc_,oo ||/ifc ||z,°°(f2) = 0, then + must be empty, otherwise 
there would exist T G T + and for all k we would have ||/ifc||i,o°(n) > IT] 1 /**. □ 

This lemma, as Lemma 3, is just a geometric observation, and a consequence of the fact that the 
sequence of triangulations is shape regular and obtained by refinement, but it is independent of the 
particular problem being considered. 

As an immediate consequence of Definition 5 and Lemma 4.1 in [MSV07] we have that 

T + = U V- 

fc>0 

Remark 3. Theorems 4.2 and 4.3 hold independently of the marking strategy. In the next theorem, we 
will make use for the first time of the assumption (2.17) done on the module MARK. 

Theorem 4.4 (Estimator's convergence: Third part). If {Tk}k&^ denote the triangulations correspond- 
ing to the convergent subsequence of discrete eigenpairs from Theorem 3.1, then the contribution of T k + 
to the estimator vanishes in the limit, i.e., 

lim m (T+) = 0. 

k — >OG 

Proof. Let T G T + , then there exists ko such that T G Tf., for all k > ko- Taking into account 
that all marked elements are at least refined once, we have that T ^ M.k- From assumption (2.17), 
m( T ) ^ Vk(Mk). Since M k C T£ U T^°, Theorems 4.2 and 4.3 imply that 

Vk(T) 2 < V k(Mk) 2 < Vk(%*) 2 + Vk(T k ) 2 — 0. 

We have thus proved that 

T)k(T) — >0, foraUTeT+. 
Now, we will prove that, moreover, 

]T r lk (T) 2 ^0. 

T6T+ 

To prove this, we resort to a generalized majorized convergence theorem. Wc first define 
e k\r ■= ypJ)k{T) 2 , for all T 6 7^", and e k := 0, otherwise. 

Then X^TeT + Vk(T) 2 = / fi ej(i) dx, and e k (x) — > as k — > oo for almost every x G f2. It remains to 
prove that J fi e k {x) dx — > as k — > oo. 

Let k be fixed. Due to the definition of T k + , for T G T k + we have that uik(T) = u)j(T) for all j > k, 
and we can drop the subscript and call this set u>(T). Using Corollary 1 we have that if (A, it) is any 



14 



fixed eigenpair of (2.1), 

m{T? < ||V(u fc - u)\\l {T) + \\\u\\l {T) + (1 + \ k ) 2 \\u k \\ 2 m(u{T]] 

< ||V(«fc - Woo)||'(T) + l|VUoo||*(D + W\l(T) + \\M\l{T) 
+ (1 + A ) 2 |K. - Uoo||li-1( W (T)) + (1 + Ao) 2 ||Uoo||fl-iMT)) 



< 



(1 + A ) 2 (JK - Moo||hi (w ( T )) + 4) > 



where 
and fulfills 



c 2 



t — ll u oo|| H i (w(T)) + ||Am|L (t) + W^uW^y 
J2 ll u oo|| 2 ffl( ^) + ||A U || 2 2 + ||V U || 2 < oo. (4.1) 



TGT+ 

Let now M k be defined by 
C 

M k\r ■= jyj (\\u k - Woo||lfi( w (T)) + c t)> for all TeT/, and M fe := 0, otherwise, 

where C is chosen so that < e k (x) < M k (x), for all x g CI. If we define 

c 2 

M\ T := C-!j-, for all TeT+, and A/ := 0, otherwise, 



then 



/ \M k (x) - M{x)\ dx = E / \M k (x)-M(x)\ dx+ E / l-^fcW ~ M ( x )\ dx 

E / i^toi E ii 

— , , l J T _ 1 



rer+\r+ ' 

<C E 4 + CIK - Woo|| 2 ^ l(n)• 
TGT+\r+ 

The terms in the right hand side tend to zero when k tends to infinity, due to (4.1) and the fact that u k 
converges to itoo in Hq(SI). Therefore, 

M fe — ► M, in L 1 (fi+). 

Hence, using that e k (x) — > 0, for almost every x £ fi, we can apply a generalized majorized convergence 
theorem [Zei90, p. 1015] to conclude that 



i\ + ) 2 = E ^m 2 = / dx — > °< 



?7fc 

as fc — » oo. □ 

We have proved in this section that r) k (T k ) — > as fc — > oo. In the next section we will use this result 
to conclude that (Aoojitoo) is an eigenpair of the continuous problem (2.1). 



5 The limiting pair is an eigenpair 

In this section we will prove that (Aooj^oo) is an eigenpair of the continuous problem (2.1). The idea 
in [MSV07] to prove that u^o is the exact solution to the continuous problem, consisted in using the 
reliability of the a posteriori error estimators, that is, the fact that the error in energy norm is bounded 
(up to a constant) by the global error estimator. Such a bound does not hold in this case unless the 
underlying triangulation is sufficiently fine (see Theorem 2.1). We do not enforce such a condition on 
the initial triangulation 7o, since the term sufficiently fine is not easily quantifiable. Instead we resort 
to another idea, we will bound a(uoc,v) — A 00 6(u 00 ,w) by the residuals of the discrete problems, which 
are in turn bounded by the estimators, and were proved to converge to zero in the previous section. 
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Theorem 5.1. The limiting pair (AqcMoo) of Theorem 3.1 is an eigenpair of the continuous prob- 
lem (2.1). That is, 

a(uoo,v) = Xoo 6(1/00,^), V v S i?o(fi), 



||Uc 



1. 



Proof. We know that ||itoo||f, = 1 due to Theorem 3.1. It remains to prove that 

a(woo, v) = Aoofrf/Uoo, v), Vue Hi (f2). 
Let v G Hq(£1), and let v k S Vfc be the Scott-Zhang intcrpolant [SZ90],[SZ92] of v, which satisfies 

1 /2 

ll u - «fc|| T < /iT||Vu|| Wfc(T ) and \\v - v k \\ dT < || Vu|j Wfc(T) . 
From (2.16) we have 



a(u k ,v k ) = X k b(u k ,v k ), 



for all k, and then 



(5.1) 



\a(u 00 ,v)-X QO b(u 00 ,v)\ = |a(itoo,^) - A 00 fe(w 00 ,w) - a(u k ,v k ) + X k b(u k ,v k )\ 

= \a(u k ,v- v k ) - X k b(u k , v - + fe(A fe it fc - XooUoo,v) + a(itoo - Ufc,«)| 
< |a(it fc ,w - u fe ) - X k b(u k ,v - v k )\ + \b(X k u k - A^Uoo, u)| + |o(itoo - Uk,v)\. 

The second term in (5.1) can be bounded as 

\b(X k u k - AooUoo, v)\ = \X k b{u k - Moo, u) + (Aft - A 00 )6(w 00 , u)| 

< \X k \\b(u k - Uoo,v)\ + \X k - A o||6(u 00 ,'u)| 

< A - woolly ||u|| n + - A oo| ||uoo|ln IMIn 

< (A - -Woolly + |A fe - Aoo| ||uoo|lo) ||«|| n ■ 

And the third term in (5.1) is bounded by 

[a(«oo - Uk,v)\ < ||V(itoo - Ufc)||n ||Vv|| n . 

Finally, the first term in (5.1) can be bounded following the steps of the proof of the a posteriori 
upper bound, as follows: 



\a(u k ,v- v k ) - X k b(u kl v - v k )\ 



/ AVu k ■ V(v - v k ) - X k / Bu k (v-v, 
V / (-V-(AVu k )-X k Bu k )(v-v k ) + / (v-v k )AVu k 

n eTk JT JdT 

[ R(Xk,Uk){v - v k ) + \ I {v - v k )J(u k ) , 

n^r- JT 1 JdT 



with R(X k , u k ) and J(u k ) as defined in (2.4) and (2.5). Now, by Holder and Cauchy-Schwarz inequalities 
we obtain 

\a(u k ,v - v k ) - X k b(u k ,v - v k )\ < \\R(X k ,Ufc)|| T \\v v k \\ T 

dT \\ v v k WdT 

TeT k 

< \\R(Xk,u k )\\ T h T \\Vv\\ u> k (T) + \\J(uk)\\g T h^WVvW^^T) 

T£T k 

1/2 



<(j2 411^,^)11^ + ^11^)11^ nv«|| 

\TeT fc / 



Vk(%) \\^v\\ n . 
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Summarizing, we have that 
|a(«oo,") - A 00 fe(u 00 ,w)| < ((1 + A ) \\uk - Wcx)||jji(n) + - Aoo| ll^oolln + Vk(T k )j IMI.ffi(n) ■ 

Using the convergence of Uk to Uoo in and Xk to Aoo in K from Theorem 3.1, and the convergence 

of the global estimator to zero from Theorem 4.1, we conclude that 

\a(uac,v) - A 00 6(u 00 , v)\ = 0, 

and the proof is completed. □ 



6 Main result and concluding remarks 

We conclude this article by stating and proving our main result, which is a consequence of the results in 
the previous sections, and discussing its strengths and weaknesses. 

Theorem 6.1. Let {(Xk, Uk)}keN denote the whole sequence of discrete eigenpairs obtained through the 
adaptive loop stated in Section 2.4- Then, there exists an eigenvalue X of the continuous problem (2.1) 
such that 

lim Afe = A and lim dist^i(Q) (uk, M(X)) = 0. 

Proof. By Theorem 3.1, taking A := A^, we have that lim^oo Xk = A, and by Theorem 5.1, A is an 
eigenvalue of the continuous problem (2.1). In order to prove that lim dist^i(m (uk, M(A)) = we argue 

by contradiction. If the result were not true, then there would exist a number e > and a subsequence 
{wfe m }meN 01 {wfcjfceNo such that 

dist H i (n) (u fcmJ M(A))>e, VmeN . (6.1) 

By Remark 2 it is possible to extract a subsequence of {uk m }meN which still converges to some function 
Uoo £ Vqo- By the arguments of Sections 4 and 5, Uoq is an eigenfunction of the continuous prob- 
lem (2.1) corresponding to the same eigenvalue A. That is, a subsequence of {wfe m }meN converges to an 
eigenfunction in M(A), this contradicts (6.1) and completes the proof. □ 

Remark 4. We have proved that the discrete eigenvalues converge to an eigenvalue of the continuous 
problem, and the discrete eigenfunctions converge to the set of the corresponding continuous eigenfunc- 
tions, and this is the main result of this article. But there is still an open question: If Afc was chosen 
as the j-th eigenvalue of the discrete problem over 7fc, is it true that {Afc}fc g N converges to the j-th 
eigenvalue of the continuous problem? The answer is affirmative for a large number of problems, but not 
necessarily for all. There could be some pathological cases in which looking for the j-th eigenvalue we 
converge to one that is larger. 

We now state an assumption on problem (2.1) that we will prove to be sufficient to guarantee that the 
convergence holds to the desired eigcnvaluc/cigcnfunction. More precise sufficient conditions on problem 
data A and B to guarantee that this assumption holds will be stated below. 

Assumption 1 (Non-Degeneracy Assumption). We will say that problem (2.1) satisfies the Non- 
Degeneracy Assumption if whenever u is an eigenfunction of (2.1), there is no nonempty open subset O 
of ft such that u\ Q G Vt{0). 

Theorem 6.2. Let us suppose that the continuous problem (2.1) satisfies the Non-Degeneracy Assump- 
tion 1, and let {(Xk,Uk)}keN denote the whole sequence of discrete eigenpairs obtained through the 
adaptive loop stated in Section 2.4 and X denote the j-th eigenvalue of the continuous problem (2.1). 
Then, 

lim Afc = A and lim dist ff i/Q) (uk, M(A)) = 0. 

k— >oo k— >oo 

Before embarking into the proof of this theorem, it is worth mentioning that the model case of A = / 
and B = 1 satisfies Assumption 1, due to the fact that the eigenfunctions of the laplacian are analytic. A 
weaker assumption on the coefficients A and B that guarantee non-degeneracy of the problem are given 
in the following 
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Lemma 6. If A is continuous, and piecewise V\, and B is piecewise constant, then problem (2.1) satisfies 
the N on- Degeneracy Assumption 1. 

Proof. We will argue by contradiction. Let us suppose that there exists an eigcnfunction u of (2.1) with 
corresponding eigenvalue A, and a nonempty open subset O of Vl such that u\o € Ve(0). Without loss 
of generality, we may assume that A\o € Vi(0) and B is constant over O. Then 

-V • (AVu) = XBu, in O. 

Since u\q G Ve(0), we have that —V • (AVu) G Ve-\{0), and the last equation implies that u\q G 
Ve-i{0). Repeating this argument we finally obtain that 

u\o = 0, 

which cannot be true. In fact, u is a solution of a linear elliptic equation of second order with uniformly 
elliptic and Lipschitz leading coefficients and therefore, it cannot vanish in an open subset of f2 unless it 
vanishes over Q, [Han94]. □ 

Remark 5. Searching for other sufficient conditions on the coefficients to guarantee Assumption 1 is 
out of the scope of this article. We believe that in the assumptions of the previous lemma, A can be 
allowed to be piecewise continuous with discontinuities along Lipschitz interfaces. The only thing needed 
is a proof of the fact that solutions to elliptic problems with coefficients like these cannot vanish in an 
open subset of ft unless they vanish over all O. We conjecture that this could be proved using Han's 
result [Han94] in combination with Hopf's lemma [GT83], but it will be subject of future work. 

We now proceed to prove Theorem 6.2, which will be a consequence of the following lemma. 

Lemma 7. Let {hk}keN denote the sequence of meshsize functions obtained through the adaptive loop 
stated in Section 2-4- If the continuous problem (2.1) satisfies the Non-Degeneracy Assumption 1, then 
HMU°°(n) -* as k -> oo. 

Proof. We argue by contradiction. By Lemma 5, if ||/ifc||z,=°(n) does not tend to zero, then Q + is not 
empty, and then there exists T G 7~ + , and thus fco G No such that T G Tk, for all k > fco- Since 
||itfc m — Uoo\\l 2 (t) - * as m — > oo, and Uk\r £ Ve(T), for all k > 0, using that Vi(T) is a finite 
dimensional space we conclude that 

u 00 \ T eV l {T). (6.2) 
Theorem 5.1 claims that Uoo is an eigenfunction of (2.1) and thus (6.2) contradicts Assumption 1. □ 

Remark 6. It is important to notice that the convergence of hk to zero is not an assumption, but a 
consequence of the fact that a subsequence is converging to an eigcnfunction and the Non-Degeneracy 
Assumption 1. 

Proof of Theorem 6. 2. In view of Theorem 6. 1 it remains to prove that A& converges to the j-th eigenvalue 
of (2.1). By Lemma 7 the result follows from (2.3). □ 

We conclude the article with several remarks. 

Remark 7. At first sight, the convergence of ||ftfc||£«>(n) to zero looks like a very strong statement, 
especially in the context of adaptivity. But the uniform convergence of the meshsize to zero should not 
be confused with quasi-uniformity of the sequence of triangulations {Tk}k£n , the latter is not necessary 
for the former to hold. Thinking about this more carefully, we realize that if we wish to have (optimal) 
convergence of finite element functions to some given function in H l (£l), then hk must tend to zero 
everywhere (pointwise) unless the objective function is itself a polynomial of degree < I in an open 
region of Q. Lemma 3 implies that the convergence of hk to zero is also uniform, and this does not 
necessarily destroy optimality [CKNS07, Ste07, GMZ08]. 

Remark 8. A sufficient condition to guarantee that we converge to the desired eigenvalue is to assume 
that hk — > as k — > oo. This condition is weaker than the Non-Degeneracy Assumption, but it is in 
general impossible to prove a priori. 
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Remark 9. Another option to guarantee convergence to the desired eigenvalue is to start with a mesh 
which is sufficiently fine. In view of the minimum-maximum principles, it is sufficient to start with a 
triangulation Tq that is sufficiently fine to guarantee that Xj,T < Xj , where jo > j is the minimum index 
such that Aj Q > Xj . This condition is verifiable a posteriori if we have a method to compute eigenvalues 
approximating from below. Some ideas in this direction are presented in [AD03], where the effect of 
mass lumping on the computation of discrete eigenvalues is studied. 
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